Density expansion for transport coefficients: Long-wavelength versus Fermi surface 

nonanalyticities 
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The expansion of the conductivity in 2-d quantum Lorentz models in terms of the scatterer density 
n is considered. We show that nonanalyticities in the density expansion due to scattering processes 
with small and large momentum transfers, respectively, have different functional forms. Some of the 
latter are not logarithmic, but rather of power-law nature, in sharp contrast to the 3-d case. In a 2-d 
model with point-like scatterers we find that the leading nonanalytic correction to the Boltzmann 
conductivity, apart from the frequency dependent weak- localization term, is of order n'^^^. 



PACS numbers; 51.10-|-y, 05.60+w 

It is well known from the statistical mechanics of fluids 
that for transport coefficients, as opposed to thermody- 
namic quantities, no virial or density expansion exists [Q. 
Let us consider a classical Lorentz gas i.e. a single 
particle moving in a static array of scatterers with scat- 
terer density n, as a simple model of a classical fluid. In 
such a system in three dimensions (d = 3), the diffusion 
coefficient, D, as a function of n has the form, 

D/Db = 1 + Din + D2inn'^\nn + Di^^ +o(n^) . (la) 

Here Db denotes the Boltzmann diffusivity, Di, D2in, 
and D2 are numbers, and o(n^) denotes terms that van- 
ish faster than for n — > 0. In 2-d systems, a similar 
nonanalyticity appears, but at one lower order in the 
density expansion. 



D/Db = 1 + DiinnlTLn + Din + o{n) 



(lb) 



The nonanalyticities in Eqs. (1) are not specific to the 
Lorentz models, but are also present in real fluids 
They are a result of long range dynamical correlations 
in the system. If one performs a cluster expansion in a 
d-dimensional system, then ring collisions, i.e. processes 
where the scattered particle collides with a scatterer to 
which it returns after having scattered off a number of 
other scattering centers, lead to a logarithmic infinity in 
the density expansion at order n'^~^ This divergence 
is cut off by the mean-free path. Through the density de- 
pendence of the latter, this translates into a logarithmic 
nonanalyticity at that order. The nature of the higher 
order terms is not known, but they are believed to also 
contain logarithmic nonanalyticities. 

Since the mechanism that is believed to lead to these 
effects in classical systems is rather general, one would 
expect a similar effect to occur in the transport coeffi- 
cients of quantum mechanical particles. This is indeed 
the case Although the leading expansion parame- 
ter is different and performing the classical limit to 
make contact with the classical ring collisions is nontriv- 
ial [Q, one obtains for a 2>-d quantum system again an 



expansion of the form given in Eq. ([la|). Specifically, all 
of the coefficients in Eq. (|l^) have been calculated for a 
quantum Lorentz model that is a good representation of 
electrons injected into Helium gas Q . It is apparent from 
the details of these calculations, although not from the 
result, that in a quantum system there are two physically 
distinct sources of the logarithmic nonanalyticity. These 
are, (1) long-wavelength contributions, i.e. those domi- 
nated by scattering processes with a momentum transfer 
on the order of an inverse mean-free path, and (2) Fermi 
surface contributions, which are dominated by momen- 
tum transfers close to 2fcF, with A;f being the Fermi mo- 
mentum. In a dilute system, these two length scales are 
well separated. Again, the nature of the higher terms 
is not known, but inspection of some individual terms 
makes it appear likely that the next term in Eq. ( p^ ) is 
of the form ri'^(lnn)^. 

In this Letter we show that the fact that both of the 
mechanisms mentioned above yield a logarithmic nonan- 
alyticity is characteristic of d = 3, and that in d = 2 
some of the Fermi surface contributions lead to a power 
law nonanalyticity of the form n'^/^. Moreover, we find 
that in a 2-d quantum Lorentz model with point-like 
scatterers, the leading nlnn nonanalyticities cancel, so 
that the n^^"^ term is the leading nonanalyticity, apart 
from the frequency dependent weak-localization loga- 
rithm that appears in a 2-d quantum system. The den- 
sity expansion for the frequency dependent conductivity, 
(which in a quantum system is easier to calculate 
than the diffusivity) in such a model thus takes the form. 
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where we have left out contributions that vanish as 
uj 0. Here (TB=e^er/(27rTO) denotes the Boltzmann 
conductivity, r is the scattering mean-free time, m is the 
electron mass, e = /cp, and 7 = m/r = nurkp, with or 
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the total s-wave cross-section, is a convenient expansion 
parameter that is proportional to the impurity density. 
The In(wr) term is the so-called weak-localization cor- 
rection (see below), whose prefactor has been known for 
some time. For the remaining coefficients in Eq. we 
find, 
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= 3/2V2. (2b) 



Equations (g) constitute our result. Before we sketch 
its derivation, let us explain the physical origin of the 
7'^/^ term in Eq. (pa|). In a degenerate system of nonin- 
teracting fermions, the Pauli principle restricts the phase 
space that is available in scattering processes. These re- 
strictions lead, e.g., to the well-known nonanalyticity of 
the Lindhard function at a momentum q = 2kp. The na- 
ture of this nonanalyticity is dimensionality dependent; 
it is logarithmic in rf = 3, and a square root in = 2. 
The same phase space restrictions lead to related non- 
analyticities in the scattering cross section, and hence 
in the transport coefficients. This is precisely what we 
find. The effect we predict is thus a consequence of the 
sharpness of the Fermi surface. It is related to other phe- 
nomena resulting from the degenerate nature of a Fermi 
gas at r = 0, like, e.g., the Friedel oscillations, and the 
Kohn anomaly. 

We now outline the derivation of our results. Let us 
consider the standard Edwards model of noninteracting 
electrons in d = 2 in an environment of static, spatially 
random scatterers. The model Hamiltonian reads. 



= ^ (ek - m) flkflk 

k 



Ev(q) 

k,q 



'k+q/2 °^k-q/2 , (3) 



where and ak denote the creation and annihilation 
operators, respectively, for electrons with wave vector k, 
fi is the chemical potential, V^(q) is the Fourier trans- 
form of the electron-impurity scattering potential, and 
Ek = k^/2TO. Throughout this paper we use units such 
that h = 1. Also, we will be working at zero tempera- 
ture, so we put ^ — ep = fcp/2m. Since we are dealing 
with noninteracting electrons, spin just leads to trivial 
factors of two and can be omitted. Standard diagram- 
matic perturbation theory j^] is formulated in terms of 
retarded (R) and advanced (A) zero temperature Green's 
functions, 
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Lj- H±iO 



(4a) 



and their impurity averaged counterparts {G^'p{ijj)}dis, 
where {. . .}dis denotes the average over the quenched 
disorder. It is most convenient to keep only the low- 
est order contribution to the self energy in the averaged 
Green's function, so we will use as the building blocks of 
our perturbation theory the following approximation to 
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(ck - ep) ± n/2m 



(4b) 



In wr iting Eq. (||), we have specialized to the case of 
pure s-wave scattering, or a point-like impurity potential, 
i.e. we have put V"(q) = U ^ -^j^-K^rr? jl^. 

The transport coefficient of interest to us, viz. the dy- 
namical conductivity cr(a;), can be expressed in terms of 
the Green's functions by means of the Kubo- Greenwood 
formula, 

Rea(c.) = ^Re ^ fc, {QlM^i^^ = 0) 
k,p 

-e^pH5^;k(^ = o)},j, p. , (5) 

This expression can be used to systematically expand tr 
in powers of 7. Such an expansion has been set up in 
Rcf . , and in Ref . ||] all diagrams were identified that 
contribute up to and including order 7^ in d = 3. The 
classification of the diagrams with respect to the order in 
7 they contribute to does not carry over to other values 
of d, but it turns out that the diagrams that contribute 
to the terms shown in Eq. (|^ form a subset of those 
considered in Ref. They are shown in Figs. 0- H. 

All of these diagrams can be expressed in terms of inte- 
grals over combinations of two functions, J"'"'" and J"' , 
that are convolutions of and . At zero frequency 
they are defined as, 



J-(.)^/dk-^ 



I 



+ 17 e — (k — q)^ + iv^ 
with = ±. Doing the integrals yields, in d = 2, 
27r 



, (6a) 



J++(g) 



[ln(w++ + g) - ln(w++ - q)\ , (6b) 
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(6c) 



where, 



W++ = \/ w{q) - 4i7 , w+_ = \/-(^w{q) + 47^ , 

(6d) 

with w{q) = — 4e. 

The only diagrams that contribute to the 7^/^ term in 
Eq. ( [2a| ) are (e), (f), and (g) in Fig. In order to demon- 
strate how the nonanalyticity arises, let us consider dia- 
gram (g) as an example. After simple manipulations, its 
contribution to the conductivity can be written. 
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FIG. 1. Simple diagrams that contribute to the terms 
shown in Eq. (^^. For each of the diagrams (f) and (j) there 
is one equivalent symmetric diagram (not shown) that also 
contributes, and the complex conjugates of diagrams (b) - (f) 
and (h) - (j) contribute as well. The 'triangulated' line in 
diagram (b) denotes a Green's function, Eq. (^b|), with its 
value at 7 = subtracted to avoid double counting. Dia- 
gram (a) yields as + 0(7^) and serves to normalize all other 
contributions. Diagrams (e) - (g) contribute to the j^^^ term. 



+ c. c. 







+ c. c. 



FIG. 2. Infinite resummations that contribute to the terms 
shown in Eq. (|2^). The complex conjugates of diagrams (b) - 
(d) contribute as well. 
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FIG. 3. Definition of the vertices in Fig. ^. 

TABLE I. Contributions of diagrams as shown in Figs. ^ 
and ^ to the coefficients in Eq. (^^. No entry means that 
the diagram does not contribute for structural reasons, while 
a zero indicates no contribution due to cancellations. 



diagram 



0"lln 



0"1 



1 (b) 

1 (c) 

1 (d) 

1 (e) 
1 (f) 

1 (g) 
1 (h) 

1 (i) 

1 Ci) 

2 (a) 
(b) + (c)- 



-2/7r 



2/7r 


-(8/7r)ln2 



Q'3/2 



-\/2/2 

V2 
V2/4 






-(d) 



2/7r 



(4/7r)ln2 
(2/7r)ln2 



Here we have put lo — 0, and have kept only the lead- 
ing contribution for 7-^-0. It is now easy to see how 
the nonanalyticity arises. Equation (6b) shows that, in 
the limit 7-^-0, contains a singularity of the form 
{q ~ 2fcF)~^/^. A g-integration over (J++)'^ thus leads 
to a 7~i/^ term. Since ( J"'""'")'^ first appears in the inte- 
grands at order 7^, the leading singularity produced by 
this mechanism is of the form 7'^/^. Asymptotic analy- 
sis yields the prefactor of the nonanalyticity, which gives 



a contribution to the number cr3/2 in Eq. (2a) as stated 
in Table |. Similarly, an integral over ( J+^^ produces 
an In 7, and this mechanism contributes to the prefac- 
tor CTi In in Eq. ( pa| ) . Another contribution to a\ in comes 
from the 'long-wavelength' terms, which manifest them- 
selves as integrals / dq \J^~{q^. Such integrals arise 
from Fig. |l| (d) (in (c) they cancel), and Fig. g (a). The 
fact that in c? = 2 two powers of J"' are sufficient to 
produce a logarithm accounts for the fact that this non- 
analyticity appears already at linear order in 7. Finally, 
some of the diagrams discussed so far, and all of the re- 
maining ones, contribute to the analytic term at order 7. 
The various contributions are listed in Table |. 

The infinite resummation denoted by diagram (a) in 
Fig. H plays a special role in our perturbation theory, and 
deserves some discussion. This 'crossed ladder' resumma- 
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tion is the only diagram where the zero-frequency hmit 
needs to be handled with some care, since it is related to 
the so-called weak-localization anomaly, i.e. the fact that 
the conductivity of disordered noninteracting electrons in 
d = 2 contains a Inwr in perturbation theory, and that 
the true zero- frequency value of a is zero Taken 
at face value, the diagram is finite at cj = 0, since the 
resummation leads to a structure 1/(1 — 7J"' (g)/7r^). 
Expanding in powers of 7 leads to a diffusion pole, 
i.e. a 1/g^ singularity, at lowest order, but the subleading 
contribution of 0(1) to J"' {q = 0) seems to protect the 
singularity. This is misleading, however, since it is well 
known that the 'crossed ladder' is an approximation to 
an exact vertex part A(q,uj) that has an exact diffusion 
pole, A(g,aj) ^ 1/ {~iuj + q'^ / D), with D the diffusion co- 
efficient . Indeed one can show that there exist classes 
of diagrams that cancel the mass in the simple crossed 
ladder, order by order in perturbation theory Al- 
though formally of higher order in 7, these contributions 
lead to the logarithmic singularity stemming from A be- 
ing protected only by a finite frequency, and not by 7. On 
a more formal level, strict perturbation theory in pow- 
ers of 7 violates a Ward identity that reflects particle 
number conservation in the presence of time reversal in- 
variance. By choosing a self energy that is related to the 
crossed ladder vertex correction by means of this Ward 
identity, one can construct a conserving approximation 
for A which has the diffusion pole built in. Replacing 
the small-g part of the crossed ladder by the appropri- 
ate exact diffusion pole then leads to the term In cot with 
a prefactor as first reported by Gorkov et al., Ref. [ pT| , 
and made rigorous by Kirkpatrick and Dorfman and 
shown in Eq. (|2^). Apart from this, the first term in the 
infinite crossed ladder resummation also contributes to 
the coefficients ai and crun, see Table |. 

We have ascertained that no other diagrams contribute 
to the terms we are considering. For simple diagrams, it 
is easy to show, by a combination of the diagram rules 
with power counting, that diagrams with more than three 
impurity lines cannot contribute. For the infinite rcsum- 
mations, the same type of argument shows that dressing 
the diagrams shown in Fig. |^ by additional impurity lines 
leads to terms that vanish faster than j^^^. Finally, ar- 
guments analogous to those employed before in d = 3 
1^ show that one need not consider diagrams with more 
than one ladder or crossed ladder resummation. 

We also mention that the temperature dependent con- 
ductivity is easily obtained from our T = result by 
a convolution with the derivative of a Fermi function 
|p^,p[. At low temperatures, the 7'^/^ anomaly then gets 
translated into a T"^/^ dependence at fixed scatterer den- 
sity. After subtracting out the weak-localization term, 
this should be observable, at least in principle, in exper- 
iments of the type reported by Adams |l^ . 
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